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Abstract: I reconsider the problem of the Newtonian limit in nonlin- 
ear gravity models in the light of recently proposed models C grav ~ 
y/—gf{R) with inverse powers of R. 

Expansion around a maximally symmetric local background with cur- 
vature scalar Rq > gives the correct Newtonian limit on length 
scales <C R if the gravitational Lagrangian y/—gf(R) satisfies 
\f(R )f"(R )\ < 1, and I propose two models with f"(R ) = 0. 
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The need for an effective or genuine cosmological constant to explain the faster than 
expected cosmological expansion in our epoch has generated a lot of activity on scalar 
field (" quintessence" ) models, where the potential energy or an otherwise anomalous 
dispersion relation of the quintessence accelerates the expansion. 
On the other hand, it is known that curvature terms can also accelerate the expansion 
of the universe [3 13 13 13 13 E| > and the application of this mechanism to explain the 
current expansion rate has been denoted as curvature quintessence jjj. While this 
mechanism usually relies on higher order curvature terms, it has also been noticed 
recently that inclusion of an term in the gravitational Lagrangian would yield a 
scale factor a(t) cx t 2 01313. 

The model proposed recently by Carroll et al. (CDTT), C ~ R— (/i 4 /i?) 0, fits into 
the framework of the so-called nonlinear gravity (NLG) models 

M 2 

£ — —^-\^gf(R) + ^-matter, (1) 

see (3 [ni] and references there, and for brevity I denote models with /(0) = oo as 
singular NLG models in the sequel. 

The generalized Einstein equations following from are 
G, u = f{R)R flv - l -f{R)g,u - V„V„/'(i2) + g^V 2 f'(R) 



and it is readily verified that V^G^ = 0. NLG theories usually assume 1 f(R) = 
R + Q£ 2 R 2 + 0(R 3 ), whence © admits flat Minkowski space as a maximally symmet- 
ric vacuum solution, and the Newtonian limit proceeds as in Einstein gravity, with 
additional Yukawa terms in the gravitational potential [121 EH EEH UH] ■ Suppression 
of the Yukawa terms at macroscopic distances leaves only the the leading 1/r term, 
and one finds that M = Mpi = (87rGjv)~ 1 / 2 is the reduced Planck mass as in Einstein 
gravity. 

However, the model proposed in [5] does not allow flat Minkowski space as a solution, 
and the problem of the Newtonian limit is more intricate 2 . Intuitively one would ex- 
pect that on length scales much smaller than an intrinsic curvature scale one should 
be able to recover the Newtonian limit, but intuition can be deceiving, and it is known 
in the framework of regular NLG models that these models may not admit a consis- 
tent weak field approximation. Therefore I propose the following approach to study 
this problem for singular NLG theories: Since our four- dimensional spacetime locally 
admits a ten-dimensional group of symmetry transformations, the Newtonian limit, if 
it exists, should be recoverable from expansion around a maximally symmetric local 
background geometry, which contrary to the regular case now will have to correspond 
to a curvature scalar R ^ 0. This will be used in Sec. |3]to demonstrate that existence 
of a weak field approximation around a symmetric local background with Ricci scalar 
Rq > can be achieved by imposing the condition f"(Ro) = on the singular NLG 
models. In these models M is then related to the reduced Planck mass through 

M = Mpi/ y/fW ■ 

However, before entering the discussion of the Newtonian limit in singular NLG mod- 
els, I would like to revisit and slightly extend the evidence for accelerated expansion 
in these models in Sec. 121 

2 The cosmological behavior at late times 

The cosmological evolution equations from (J2J) are quite complicated, but we can 
make a general statement about the late time expansion behavior of singular NLG 
models. 

1 I follow the MTW conventions ^I] f° r the signature of the metric and the definition of the Ricci 
tensor: 

p per f) f) per _i_ per pp p(T pp 

r^pv — /iffy — Ufjl. fji* U V L ^(j T" -L per 1 pv 1 pv^ pa- 
It is useful to keep that in mind when comparing with the literature on regular NLG models, because 
the relative signs between even and odd powers of R depend on these conventions. 

2 Capozziello et at had noticed that f(R) — R^ 1 yields a(t) oc t 2 0|E], but did not further pursue 
this model. f(R) — R^ 1 would not have a maximally symmetric vacuum solution. 



Since the generalized Einstein equation (J2J) still implies energy-momentum conserva- 
tion V^T^, = 0, the time evolution of the scale factor a(t) in a Friedmann model is 
still governed by the generalized Friedmann equation SC/5g 00 \FRw metric = 0. 
For the spatially flat FRW metric 

ds 2 = -dt 2 + a 2 (t)dx 2 

the generalized Friedmann equation following from (j2J) is 

„. / d d 2 \ d 1 „ / d d 2 \ d _ „, / a d 2 
-3f 6- + 6- - + -/ 6- + 6- + 3-dof 6- + 6- 
\ a cr J a 2 \ a cr J a \ a er 

with T 00 = Q- 

In general this will be a third order equation for the scale factor. To analyze the late 
time behavior, we first assume that there is only ordinary dust and radiation in g, 
whence the energy density can be neglected at late times for expanding solutions. 
We then make a power law ansatz a(t) oc t a , which yields 

-3/'g a (2a-l))^ + i/g a( 2 a -l) 
_ 3( V(2a-l) / „/6 N a 



t 4 J \t 2 

If R~ n , n > 0, is the leading order singularity in the singular NLG model f(R), then 
at late times the contribution from this term will dominate all 3 terms in Eq. (J3J), 
with the same time dependence oc t 2n . This yields the algebraic equation 

a — 1 1 n(n + 1) 

n— r H = 0, 

2(2a - 1) 2 2a - 1 

which determines the expansion coefficient a: 

2n 2 + 3n + 1 

a = . 5 

n + 2 V ; 

This was found for f(R) = R — n 2n+2 R~ n in 9J through conformal transformation to a 
corresponding scalar quintessence model, and the corresponding result for f(R) = R n , 
n > 0, is also spelled out in [16j. 

Note, however, that Eq. Q is also compatible with exponential expansion a(t) oc 
exp(Ht) at late times if f(12H 2 ) = QH 2 f'(12H 2 ) has a solution. Carroll et al. found 
in the metric formulation of their model that power law expansion is dynamically 
preferred 9j. Vollick looked at the Palatini formalism in the CDTT model and con- 
cluded that exponential expansion would arise in that formulation 17 . Our use of 
a symmetric local background geometry in the next section does not predicate the 
global late time expansion, but only assumes that spacetime should have maximal 
symmetry locally. 



3 Expansions around maximally symmetric vacua 



In the spirit of the philosophy outlined in Sec. ^we now assume that the Newtonian 
limit should be recoverable through weak field expansion around a symmetric local 
background geometry: The maximally symmetric vacuum solutions satisfy 

R afJ ,f3u = —(9a/39fiu — 9av9p.fi)-, 
R — — 

-tipu ^ 9jiV1 

where the constant curvature scalar R satisfies 

f'(Ro)Ro = 2f(Ro). (6) 

In ordinary NLG theories this always permits Rq = 0, but for singular NLG models 
this yields values Rq ^ 0, and the challenge is to derive the Newtonian limit from the 
weak field expansion around the vacuum solution. 

The first order expansion of Eq. (j2J) around a symmetric vacuum solution yields 
f(R )SR IIU + ^[f"(Ro)Ro - 2f'(Ro)]g^8R 

-\f(Ro)5g,u - f"{Ro)(V,V v 8R - g,„V 2 5R) = -j^T, u (7) 

or 

f(Ro)SR^ - ^f"(R )Rog,JR - \f{R ) - \g»Jg) 

-f(R ) (vj7 v 5R + \g^ 2 5R>j = (t, u - Ig^Tj . (8) 

The first order variation of the Ricci tensor is 

SR^ = i(V M V CT 5^ + V v V a 5g^ + 5g au R a ^ + 5g afl R a u ) 

-5g a! 3R a /» - \v 2 5g^ u - iv M V„fy (9) 

= ^(V M V ff ^ + V v V a 6g^) + ^RoSg^ - ^Rog^Jg (10) 

-^V 2 %, - -V„V„fy, 

The mass term f(Ro) ~ 0(Rq) vanishes in regular NLG theories, and should also be 
negligible in singular NLG theories at least up to length scales where Newton's law 
has been verified, which implies that Rq ~ /i 2 must correspond to a small mass scale 
\i. We also note that the mass and derivative terms following from Eqs. (jSJ) and fjlUjl 



have the correct signs for non-oscillatory attractive solutions if R$ > 0, f'(Ro) > 
(=*► f(Ro) > 0), and f"(Ro) > 0. 

In regular NLG theories f{R ) = /(0) = 0, f'(R ) = 1, and f'(R ) = 12£ 2 is assumed 
to be very small, such that the corresponding Yukawa terms are suppressed relative 
to the leading 1/r term at macroscopic distances. On the other hand, every set of 
observational tests of Newton's law can only cover a finite range of length scales. 
Therefore one might be tempted to conclude that very large £ is another possibility, 
such that e.g. the Yukawa term exp(— r/£)/r from f(R) = R + 6£ 2 R 2 at observational 
distances also approximates a 1/r term and only rescales the ratio between M and 
Mpi by a constant factor. 

That this latter possibility is excluded in regular NLG theories was noticed already 
by Pechlaner and Sexl: /"(0) = 12£ 2 has to be small for consistency of the weak field 
approximation, because otherwise domination of the fourth order terms would yield 
strong curvature on all length scales [P2] . 

This reasoning carries over to the singular case, with minor modification: Due to the 
presence of a small mass term the Newtonian potential, if it exists in the theory, will 
always come from a limit of Yukawa terms. Yet we still have to confine the impact from 
the fourth order terms to small r. This will constrain the parameter space, because 
in singular NLG theories f'(Ro) ~ and f"(Ro) ~ fi~ 2 generically would imply 

that the fourth order derivative terms dominate the equation for Sg^, thus spoiling 
the consistency of the weak field approximation. 

The need for suppression of the fourth order terms can also be seen from the following 
simple example: 

AU(r) - ^ 2 U(r) - ^A 2 f/(r) = J^tf(r) 
yields 

TYl 

U(r) = = |exp(— k^r) — exp(— k+r)] , 

87rM 2 r v /m 2 -2 / u 2 1 

with 



k± = m 2 ± my / m 2 — 2fi 2 . 
This will give a Newtonian 1/r potential at distances r only if /i m: 



U(r) w - g 7r ^ f2r exp(-/ir) - exp(-v / 2' 



mr 



1 -i -i 
, m Cr<^ . 



8ttMV 

In the terminology of the singular NLG models this means that we need 



|/(i*b)/ w (J2b)|<l, 



while e.g. f(R) = R- fi 2n+2 R- n would yield \f(R )f"(R ) \ = n(n + l) 2 /(n + 2) 2 . 
Therefore we either have to invoke a second small parameter in f(R) such that both 
f(Ro) and f'\Ro)/r 4L are small relative to f'(Ro)/r 2 at the length scales of interest. Or, 
since f(Ro) 7^ by Eq. ©, we arrange /(-R) such that the coefficient of [T 2 vanishes 
altogether, i.e. by choosing our model such that the solution of Eq. (jHJ) satisfies 

f"(Ro) = o. (li) 

In that case the fourth order terms vanish in the weak field expansion and all the cur- 
vature contributions to Eq. (JHJ) are subleading, such that for r <C fi^ 1 a flat ansatz can 
be used to determine the local potential at these scales. In leading order this is then 
nothing but the ordinary Newtonian limit at these scales, with the only modification 
that 

M = (87rG N f'(R ))-^ 2 = M Pl /vTW • (12) 

4 Two simple examples of singular NLG models 

4.1 The criterion (JTTjl is satisfied e.g. by 
„ M 2 , / R 2 3//\ „ , . 

£ = -TrV-g R + 7rT TT + Matter- (13) 



2 v a V V R 

This corresponds to 

R = 3/i 2 
and 

M = Mpi I V2 , 

and the power law for late time expansion would be the same as in the original CDTT 
model, a{t) oc t 2 . 

4.2 Another model that satisfies the criterion (jllj) is 
M 2 / u 4 u 6 \ 

£ = ^rV^g [R - 15^7 + 25 T^ + Matter- (14) 



2 v * V # 2 
This yields 

i? = 5/i 2 
and 



M = M Pi 



'5 
6' 

The R~ 2 term accelerates the power law expansion at late times to a(t) oc t 3,75 . 



5 Conclusions 



The problem of existence of a weak field expansion and the Newtonian limit is more 
intricate in singular NLG models than in regular NLG models, but can apparently 
be solved. 

Models can in particular be chosen to satisfy the constraint (|11|) to ensure consistency 

— 1/2 

of the weak field expansion at length scales <Cfl . 

Two minimal extensions of the CDTT model which satisfy this constraint are given 
in Eqs. (fT5|) and flljl . 
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